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ABSTRACT 

We report on the direct detection and characterization of the probable red supergiant progenitor of 
the intermediate-luminosity Type Il-Plateau (H-P) supernova (SN) 2012aw in the nearby (10.0 Mpc) 
spiral galaxy Messier 95 (M95; NGC 3351). We have identified the star in both Hubble Space Telescope 
images of the host galaxy, obtained 17-18 yr prior to the explosion, and near-infrared ground-based 
images, obtained 6-12 yr prior to the SN. The luminous supergiant showed evidence for substantial 
circumstellar dust, manifested as excess line-of-sight extinction. The effective total-to-selective ratio 
of extinction to the star was Ry « 4.35, which is significantly different from that of diffuse interstellar 
dust (i.e., Rv — 3.1), and the total extinction to the star was therefore, on average, Ay ~ 3.1 mag. 
We find that the observed spectral energy distribution for the progenitor star is consistent with an 
effective temperature of 3600 K (spectral type M3), and that the star therefore had a bolometric 
magnitude of —8.29. Through comparison with recent theoretical massive-star evolutionary tracks we 
can infer that the red supergiant progenitor had an initial mass 15 < Mini(M0) < 20. Interpolating 
by eye between the available tracks, we surmise that the star had initial mass ^ 17-18 Mq. The 
circumstellar dust around the progenitor must have been destroyed in the explosion, as the visual 
extinction to the SN is found to be low {Ay = 0.24 mag with Ry ~ 3.1). 

Subject headings: supernovae: general — supernovae: individual (SN 2012aw) — stars: late-type — 
stars: evolution — stars: fundamental parameters: other — galaxies: individual 
(Messier 95, NGC 3351) 



1. INTRODUCTION 

Supernovae (SNe) are among the most powerful ex- 
plosions in the Universe. In addition to Type la SNe, 
which arise from the thermonuclear runaway explosion 
of a mass-accreting white dwarf star, there are SNe 
that result from the collapse of the core at the end- 
point of a massive (with initial mass Mini ^ 8 Mq; e.g., 
IWooslev &: Weaveill 19861 ) star's evolution. If the star ex- 
plodes with most of its extended hydrogen envelope still 
relatively intact, the ev ent will be observed as a Type II - 
Plateau SN (SN II-P; IBarbon. Ciatti. fc Rosind [l97l . 
We would expect such a progenitor star to be in the red 
supergiant (RSG) phase at the time of core collapse. 

We have been extremely fortunate in recent years to 
detect and characterize the probable RSG progenitors of 
SNe II-P in nearby galaxies. (We note that the most fa- 
mous progenitor identification, of the star Sk —69° 202 
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that exploded as SN 1987A in the Large Magellanic 
Cloud, was actuall y a blue supergiant, not a RSG; e.g., 
lArnett et al.l 11989) . One of the best examples is the 
identification in high-quality, ground-based imaging data 
of the RSG progenitor of the SN II-P 2008bk in NGC 
7793 (IMattila et all \200§. IVan Dvk et all [20T2h . Other 
SN II-P RSG progenitors have also been directly identi- 
fied in archival Hubble Space Telescope (HST) images 
of nearby host galaxies , including SN 2 003gd i n M74 
(jVan Dvk. Li. fc Fihppe nko 2003; Smart t et al.l |2004), 
SN 2004A in NGC 6207 fflen drv et al. 200l), S N 2005cs 
in M 51 (jMaund. Smartt. fc Danziger 200 51: iLi et al " 



2006), and SN 2009md in NGC 3389 l|Fraser et al. 



20TH ). All five of these SNe II-P are of low lumi- 
nosity, with bolometric luminosities Lboi ^ lO^"'^'^ erg 
s~^ at maximum, and with lower ejecta velocities dur- 
ing the plateau and lower luminosity on the light-curve 
tail as a res ult of a smaller ^ ^Ni yield in the explo- 
sion (e.g., Za mpieri et al.l 120031 : .Pastore llo et al.i r2004). 
This is relative to intermediate-lumi nosity SNe H-P, such 
as SN 1999em in NGC 1637 Ce.g.. IHamuv et al.l [200ll: 
[Leonard et al.ll2002l: lEh^hamdi et al.ll2003D . with Lboi ~ 
2^q41.5_-[^q42 gj,g g-1 g^-j. niaximum. 

Three other SNe II-P also have direct probable 
progenitor identifi c ation s: SN 1999ev in NGC 4274 
Maund fc S^Mrttl [20051). SN 2004et in NGC 6946 
Li et al. 2005; Crockett et al. 2011'), and SN 2008cn 
in NGC 4608 (Elias- Rosa et al.j ,20091) . SN 1999ev has 
no published photometry or spectroscopy, so its nature 
has not been well determined. SN 2008cn appears to 
have been similar to high-luminosity SNe II-P (with 
Lboi ~ 10^^-^ erg s~^ at maximum), such as SN 1992II 
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(iClocchiatti et all Il996f ) and SN 2007od ()Inserra et alj 
1201 ID . which also exhibit a less pronounced plateau 
and more linear post-maximum decline. Furthermore, 
lElias-Rosa et alJ ([2(309) detected a pro genitor for SN 
2008cn which was more yellow than red. M aguire et aD 
(pblO. ) showed that SN 2004et bears similarities in both 
its photospheric expansion velocity and the overall shape 
of its bolometric light curve to those of SN 1999em (al- 
though SN 2004et may have been a factor of two more lu- 
minous than SN 1999em), implying that SN 2004et may 
also have been of intermediate luminosity. However, the 
nature and initial niass of the ide i itified progenitor star 
has been debated (|Li et all 120051: iCrockett et"all 120111: 
Van Dyk & Jarrett, in preparation). The progenitors of 
other intermediate-luminosity SNe II-P, therefore, have 
not yet been directly identified, although not without ef- 
forts in analyzing HST dat a to do so in the recent past, 
such as for SN 1999em (,Smartt et al.l[200l . 

The recent discovery of SN 2012aw in M95 (NGC 
3351) by P. Fagotti on 2012 March 16.86 (UT dates are 
used throughout this paper), A. Dimai on 2012 March 
16.84, and J. Skvarc on March 17.90 (reported in CBET 
3054) has now afforded us with the best opportunity 
to do so. Indications of the nature of the SN during 
the plateau phase, which we will present in a forth- 
coming paper, are that it is similar to SN 1999 e m. A 
spectrum of the SN on March 17.77 by iMunaril (|2012[ ) 
showed a very hot, blue, essentiall y featureless contin - 
uum. A spectrum on March 18.77 by Siviero 
also showed a blue, featureless continuum, and later 
spectra on March 19.85 and 19.92 exhibited the char- 
acteristics of a very young SN ITP, with the onset of 
broad lines having P-Cygni-like profiles. The object was 
also discovered b y the Palomar Tr ansient Factory (PTF; 
Eaw ct al. 2001; iRau et all l200i ) and given the name 
PTF12bvh. It was first detected at R ^ 14.2 mag in 
an im age obtained with the PTF camera (|Rahmer et al.l 
I2008f) on the Palomar Oschin Schmidt 48-in (1.2-m) tele- 
scope on March 16.70. The SN was not visible up 
to and including March 14.75 to a limiting magnitude 
i? > 22 (3(t), providing a tight constraint on the ex- 
plosion epoch (|Poznanski et al. 2012a originally reported 
R > 20.7 mag, but this has since been revised to the 
fainter limit). The SN shows indications of interac- 
tion with circumstellar matter, through detection in the 
radio (|Yadav. Chakraborti. fc Rav 2012:' Stockdale et al.l 
|2012( ) and X-ray (Immlcr & Brown 2012) bands. It also 
exhibits signs, although preli minary, of possibly un usual 
polarization of the SN ejecta ([Leonard et al.|[2012[ ). 

The probable progenitor of SN 2012aw/PTF12bvh was 
first detected in HST images by lEhas-Ro sa et al.l (|2012[ ) 
and subsequently bv iFraser et a l. (2012aj)! The appar- 
ently red color of the star indicated that it was most 
likely a RSG. Here we present photometry of the star 
and discuss its nature, including its likely initial mass. 
An analysis of the pr ogenitor has also been conduct ed by 
iFraser et al.1 ()2012bf ). In their study, iFraser et al.l con- 
clude that the RSG was observed through considerable 
visual extinction {Ay > 1.2 mag), and their estimates of 
the star's luminosity and initial mass span a fairly large 
range, lO^ ^-lO^-^ Lq and 14-26 Mq, respectively. 

For the distance to the SN, we adopt the reddening- 
and metallicity-corrected distance modul us iiq = 30.00 i b 
0.09 mag for M95 determined bv .Frecdman et al.l (|2001[) . 



2. PROGENITOR OBSERVATIONS 

The SN site was imaged by HST in F439W, F555W, 
and F814W with the Wide-Field Planetary Camera 2 
(WFPC2) between 1994 December and 1995 January, 
as part of the HST Key Project to use Cepheid-based 
distances to measure the value of the Hubble constant 
(GO-5397; PI J. Mould). The Hubble Legacy Archive 
(HLA) project has subsequently combined the indi- 
vidual WFPC2 e xposu res, usi ng the task MultiD rizzle 
(jKoekemoer et al.ll2003t see also lFruchter fc Hookl l20021. 
into deep image mosaics in F555W and F814W. The to- 
tal exposure times are 30130 s and 9830 s in these two 
bands, respectively. The HLA did not produce mosaics 
at F439W; the total exposure time is 5000 s for these 
four images. The SN site was also imaged in F555W on 
1995 December 4 by program GO-5972 (PI J. Mould) for 
2000 s total, and in F336W and F658N on 2009 January 
18 by program GO-11966 (PI M. Regan) for 4400 s and 
1800 s total, respectively. 

Images of the SN were obtained on March 21.31 in 
the i' band in sets of exposures with times of 5 s and 
10 s, under very good observing conditions 0'.'8), us- 
ing the MegaCam on the 3.6-m Canada-France-Hawaii 
Telescope (CFHT), to initially locate the SN position in 
the pre-SN archival HST images (see below) . From these 
images we also measured an absolute position for the 
SN of a(J2000) = 10h43™53!73, (5(J2000) = -Hll°40'17'.'9 
(ibO'.'ll root-mean square), relat ive to 12 stars in the field 
from the USNO Bl.O Catalog ([Monet et al.i r2003'). To 
further refine the SN position, relative to the HST im- 
ages, we obtained three sets of dithered Near-Infrared 
Camera 2 (NIRC2) images of the SN in the Kp band 
(with exposure times 1.5, 3, and 10 s, respectively), us- 
ing adaptive optics (AO) and the 10-m Keck II telescope 
on March 27.42. Since the SN was quite bright at the 
time of these observations, it was possible to perform 
the AO using it as a natural guide star, so that the laser 
guide star was not necessary. 

Using the CFHT images, we were able to locate a can- 
didate progenitor star for the SN in the HLA F814W 
image mosaic, with an error circle radius of ^ 0.6 pixel 
(O'.'OS for these drizzle d mosaics), the results of which 
were first presented bv lElias-Rosa et "ahl ([20121 ). For the 
Keck AO imaging, from a coaddition of 40 individual 3 s 
frames (no geometric distortion correction was applied 
prior to the coaddition) in which the SN is only mildly 
saturated, we used 16 stars in common between this 
coadded mosaic and the HLA F814W mosaic with the 
IRAF task geomap to register the two mosaics with an 
uncertainty in the image transformation of Aa; = 0.204 
and Ay = 0.591 HLA mosaic pixels (O'.'OlO and 0'.'030, 
respectively). We find that the location of the SN and 
progenitor candidate are consistent to within the la sta- 
tistical uncertainty of our alignment procedure. We 
conclude that we have confirmed the candidate star seen 
in the WFPC2 mosaics as the probable progenitor of SN 
2012aw. This is th e same star that was identified by 
IFraser et all ([2012al |bl). 

Using Dolph ot v2.0, as applied to WFPC2 data 
([Dolphinll2000ai rbl) , we measured the apparent brightness 
of the progenitor from the ensemble of individual im- 
ages in both F555W and F814W from GO-5397. For 
F555W we ultimately input into Dolphot only the im- 
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ages which produced uncertainties in the photometry 
of < 0.28 mag; in this band, the star is detected at 
only ~ 4(7 per individual exposure in the first place. 
We also omitted one F814W exposure which was ap- 
preciably noisier than all the others (its exposure time 
was only 230 s, compared to 1000-1500 s for the other 
exposures in this band). We used one of the remain- 
ing F814W exposures (the star is generally detected at 
22(T per exposure) as the astrometric reference image 
in Dolphot. The output from Dolphot automatically in- 
cludes the transformation from flight-system F55 5W and 
F8 14W to the corresponding Johnson- Cousins (jBesselll 
Mm) magnitudes in V and Ic- (We refer to Ic as / 
hereafter.) We find that mpsssw = 26.49 ± 0.07 and 
'TIF814W = 23.39 ± 0.02 mag. These flight-system mag- 
nitudes transform to ^ = 26.59 and / = 23.44 mag. Al- 
though our measurement at F814W is identical to that 
obtained by iFraser et al.l (j2012b.l . our measurement at 
F555W is brighter by > la than their value. For F336W 
and F439W the star was not detected by Dolphot to 
'71F336W ^ 22.4 and mF439-w ^ 25.9 mag (both 3ct). The 
measurements of the progenitor's apparent brightness (or 
limits to the apparent brightness) are summarized in Ta- 
ble [II 



As 



Fraser et al.l (|2012b[ ) discuss, there are also archival 



ground-based near-infrared images which contain the 
SN progenitor, namely those obtained at the European 
Southern Observatory (ESO), with the Infrared Spec- 
trometer and Array Camera (ISAAC) on the 8.2-m Very 
Large Telescope (VLT) Unit Telescope 1 on 2000 March 
26 and 27 (PI F. Bresolin), and with the Son of ISAAC 
(SOFI) on the 3.6-m New Technology Telescope (NTT) 
on 2006 March 24 (PI J. Ascenso). The ISAAC images 
we used were in the J, band, with 30 s individual frame 
times and four subintegrations (coadditions) in mem- 
ory, and the SOFI images were in the Ks band, with 
8 s frame times and 15 subintegrations. The progenitor 
site is also in Kg images of the host galaxy, obtained 
as commissioning data on 2000 March 23 with the Isaac 
Newton Group Red Imaging Device (INGRID) on the 
W HT. The resul t s of t hat imaging have been reported 
b vlKnaoen et all (l200l . 



j yl — _ .. . 

iFraser et al.l ( 2012b[ ) calibrated their photometric anal- 
ysis of the ISAAC and SOFI datasets using 2MASS stars 
in the images. The problem with this approach is that 
these stars are all near the photometric limit of the 
2MASS survey, J w 16.4 mag for the ISAAC images 
and Kg ~ 14.7-15.5 mag for the SOFI images. At this 
limit we expect the photometry to be far less reliable, 
as described in the 2MASS Explanatory SupplemenlQ, 
and therefore calibration of the ISAAC and SOFI im- 
age data using these stars could suffer from systematic 
effects, leading to a possible error in the final apparent 
magnitudes for the progenitor. 

We therefore reanalyzed the ISAAC and INGRID im- 
ages, including the original calibration data for those 
observations, and reanalyzed the SOFI images, apply- 
ing the calibration at Kg from the INGRID observations 
to those with SOFI, employing five well-detected stars 
in common between the two datasets. This calibration 
should be valid, since the bandpasses of both the IN- 
GRID and SOFI Kg filters are quite similar. The data 

http: / / www.ipac.caltech.edu / 2mass / releases / allsky/doc/explsup. 



reduction of the raw frames followed standard procedures 
for near-infrared imaging. For the ISAAC, SOFI, and 
INGRID data we first corrected the individual images 
for the instrumental response with a combined, normal- 
ized fiat frame obtained in each band. For the ISAAC 
observations, twilight-sky frames had been obtained for 
this purpose; for SOFI observations, dome fiats were ob- 
tained; and for the INGRID observations, we produced a 
flat from a median-filtered combination of the off-source 
sky frames. We then subtracted from the on-source im- 
ages a median-filtered combination of the sky frames. 
For each band a reference frame was chosen, all other 
frames were shifted in pixel space relative to the refer- 
ence, and the shifted frames were coadded, together with 
the reference frame, to produce a single image mosaic. 
For the March 26 ISAAC data we coadded 16 individual 
Js frames, and for the March 27 we coadded 17 frames. 
For the SOFI Kg imaging we coadded all 5 on-source 
frames. For the INGRID imaging we coadded 9 of the 
20 s frames that contained the progenitor site. 

We extracted photometry from all of the image 
mosaics using point-spread function (PSF) fitting in 
IRAF/DAOPHOT (Stetson 1987). The progenitor star 
was detected in the Jg mosaics from the two ISAAC 
nights at a signal-to-noise ratio (S/N) of ^ 25-35, 
whereas the star was detected in both the SOFI and 
INGRID mosaics at only S/N » 4. We measured aper- 
ture photometry using IRAF for the calibration-star ob- 
servations obtained for both the ISAAC and INGRID 
runs. For ISAAC these consisted of various IPersson et al.l 
(1998) standard stars observed on both nights. For both 
2000 March 26 and 27 we found solutions with airmass 
corrections of 0.100-0.105 mag airmass"^, esse ntially the 
canonical correction at J (jPersson et al.l ll998': the root- 
mean square uncertainty in the solution was 0.002 mag 
for March 28 and a slightly higher 0.009 mag for March 
27). We therefore measured J = 21.03 ± 0.04 mag on 
March 26 and J = 21.01 ± 0.03 mag on March 27 for 
the progenitor. We note that, although the images are 
in Jsi while the photometric standard magnitudes are 
in the Las Ca mpanas Observatory (LCO) Jlco band 
(jPersson et al.l[l998[ ) , we analyzed the synthetic photom- 
etry of three of these standard stars (P041-C, P177-D, 
and P330-E), for which there are also HST calibration 
spectra, and the magnitudes for these stars are only 
about 0.02 mag brighter through Js than Jlco! we also 
find that for the likely spectral energy distribution (SED) 
of the probable progenitor star (see below), the progen- 
itor itself would be about 0.05 mag brighter in Jg than 
in Jlco- (We hereafter refer to Jlco merely as J.) We 
have added this latter difference as an additional uncer- 
tainty in quadrature with the uncertainty in the photo- 
metric measurement. The uncertainty- weighted mean of 
the two measurements is J = 21.02 ± 0.03 mag. 

For the INGRID calibration, two United King- 
dom InfraRed Tel escope (UKIRT) faint standards 
(|Leggett et al.ll2006[) . FS14 and FS21, had each been ob- 
served at the low airmasses which bracketed the observa- 
tions of M95. We obtained aperture photometry of these 
two stars and calibrated the faint DAOPHOT photome- 
try of the progenitor star accordingly. We note that the 
Kg bandpass used for these observa tions is in the Mauna 
Kea Observatory (MKO) system (jSimons &: Tokunagal 
ht M)02t iTokunaga et al...2002,) . The star's brightness in 
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this band is then Ks — 19.52 ± 0.27 mag. Applying the 
cahbration of the INGRID data to the SOFI image mo- 
saic, as previously described, we obtain Kg = 19.42±0.29 
mag for the progenitor. The uncertainty-weighted mean 
of these two measurements is Kg — 19.47 ± 0.19 mag. 

We summarize the measurements of the progenitor's 
near-infrared brightness in Table [TJ Although, in the 
end, our m easurements at J a nd Kg agree with those 
obtained bv lFraser et al.l ()2012b[) to within the uncertain- 
ties in their photometry, we were able to achieve a higher 
precision in our measurements and were able to reduce 
the uncertainties substantially through the standard-star 
calibrations. 

3. THE NATURE OF THE PROGENITOR STAR 

Unlike the lower-luminosity RSGs which are progen- 
itors of lo w-luminosity SNe II- P, such a s SN 2008bk 
(iMattila et al. ^008; Van Dvk et al.ll2012[ ) , iMassev et all 
(|2005D pointed out that we would expect higher- 
luminosity RSGs to lose mass at a significantly higher 
rate, and that this mass loss results in circumstellar shells 
where dust tends to form. We already have an indication 
that this RSG was of high luminosity, from its K-h&wd 
brightness; even neglecting extinction (which, in the first 
place, at K is typically nearly a factor of 10 less than 
at V) and assuming a bolometric correction BCk ~ 3 
mag, the bolometric magnitude for the star would be 
Mboi w -7.8. 

What is striking about this progenitor detection is that 
the star is sufhciently luminous to be easily visible in 
both the HST F814W image and the ground-based near- 
infrared images. Another interesting facet is the relative 
isolation of the progenitor in the host galaxy; it is ap- 
parently far from any noticeable cluster, OB association, 
or H II region. Most likely, this isolation largely con- 
tributes to the detectability of the star at the distance 
of the host galaxy, since the star's environment is un- 
crowded. Furthermore, it is tempting to speculate that, 
like the famous Galactic RS G Betelgeuse (a Orionis; e.g., 
iNoriega-Crespo et al]|1997[ ). the SN 2012aw progenitor 
could also have been a stellar runaway. 

3.1. Variability of the Progenitor 

It is known t hat RSGs can be semiregular o r irregular 
variable stars. iKiss. Szabo. fc Beddind ()2006[ ) analyzed 
the optical variability of 48 Galactic RSGs and found 
evidence for two modes of variability in 18 of the stars 
(e.g., a Ori), one mode with periods of a few hundred 
days and one of a few thousand days. Concern therefore 
should exist whether the apparent brightnesses that we 
have measured in the various bands, especially the en- 
semble from the HST images, are representative of the 
star's actual brightness. Whereas variability at V can be 
considerable, up to a magnitude o r more, RSGs show es - 
sentially no variability at K fe.g.. iLevesaue et "al]|2007[ ). 

The HST data for M95 were obtained over a number 
of epochs spanning ^ 65 d (those data were originally 
taken to discover and measure the periods of Cepheids 
in the host galaxy), so we can analyze the multi-epoch 
photometry for any indication of short-period variability 
of the progenitor in 1994/1995, ^ 17-18 yr prior to the 
explosion. In Figure [5] we show the individual measure- 
ments with Dolphot from the ffS'T/WFPC2 images in 
both F555W and F814W. For HST observations which 



were made on essentially the same Julian date, we have 
computed an uncertainty-weighted mean value for the 
star's brightness and show those in the figure, rather than 
the two individual measurements for a given date. This 
was particularly the case for the F814W measurements. 
We also show estimates of the color from three epochs. 

Clearly excursions of the measurements exist, relative 
to the mean brightness for the star in each band, with 
total apparent variations of Ampsssw ~ 1 mag and 
AmF8i4w ~ 0-18 mag. We must determine the signifi- 
cance of these variations. For that we have computed the 
reduced statistic in each band, and find for F555W, 
;^2^^ = 1.92, and for F814W, = 5.16. We can also 
assess the p value of these two statistics; a p value of 0.05 
or less is usually regarded as statistically significant (i.e., 
the observed deviation from the mean is significant). For 
F555W we find p = 0.033, and for F814W, p = 0.001. We 
therefore consider it likely that the variab ility we observe 
is real. We note that iFraser et al.l ()2012b[ ) seemingly dis- 
miss the possibility of variability in either band. 

The mean color is already ~ 0.8 mag redder inV — I 
than the color of an early-M supergiant in the absence 
of dust. The variation in the star's color \sV — I ~ 0.6 
mag, relative to the mean. We clearly do not have spec- 
tra of the star, so we cannot determine whether the color 
variation represen ts true changes in t he star's effective 
temperature fe.g.. IMassev et alll2007f) . or merely varia- 
tions in the line-of-sight dust. If the former, the color 
variations would correspond to variations in the mean 
effective temperature of ^ 300 K. 

Although the progenitor shows considerable variabil- 
ity over the ~ 65 d at ~ which would be consis- 
tent with variations in th e dust content of the star 's cir- 
cumstehar material (e.g.. IMassev et aLll2005l [2007D . and 
more limited variability at ^ /, we will assume that the 
mean brightness at both V and / over this timescale 
is representative of the star's brightness over far longer 
timescales, at least through the year 2000, when the J 
and Kg brightness were sampled. We would expect the 
Kg brightness to have been relatively constant. We have 
no insight into the level of variability at J, but we might 
expect it to be considerably less than at /. We therefore 
assume that the progenitor's observed brightness in all of 
these bands fairly represents the actual SED of the star. 

3.2. Metallicity in the SN environment 

Assuming the oxygen abundance can serve as a proxy 
for the overall metallicity in the host galaxy, we esti- 
mated the metallicity at the SN site based on the gra- 
die nt in the disk of M95 o f the O abundance measured 
bv iMoustakas et al] (|20100 . We deprojected the Mega- 
Cam image, assuming the v alues for the posi tion angle 
and inclination for M95 from IMoustakas et al.l . and mea- 
sured the radial offset of the SN position from the nu- 
clear position. For a plate scale of 0'.'185 pixel"^, we find 
that this offset is p = 162'.'02, or 2.'70. Again, assum- 
ing the radial off s et at 25 B mag arcsec^^ for M95 from 
IMoustakas et al.l (|2010f) . i.e., p25 = 3f71, we then calcu- 
late that p/p25 = 0.73. A t this scaled nucl ear offset, as- 
suming the abundance in IMoustakas et ah] derived from 
the st rong-line index calibration bv IPilvugin fc ThuanI 
((2005h . we find that 12 + log(0/H) ?» 8.5 in the SN en- 
vironment. Given that the solar value is 12 -|- log(0/H) 
= 8.66 ± 0.05 (Asplund ct al.. .20051 we infer that the 
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metallicity at the SN site is only slightly subsolar and 
likely still consistent with solar, given the uncertainties. 
We therefore analyze our results assuming solar metal- 
licity. 

3.3. Properties of the Star 

From the relatively high apparent brightness at 
for the star, compared to the relative faintness at V, 
we realized that the reddening to the pro genitor had to 
be hig h, likely due to circumstellar dust. iMassev et aD 
(|200El) demonstrated that the effective total-to-selective 
ratio of absorption, R'y, for the dust around Galactic 
RSGs should differ considerably from the Ry = 3.1 
typical of the diffuse interstellar medium (ISM; e.g., 
ICardelli. Clavton. fc Mathi^ [Tosl . This may imply a 
difference, for example, in the grain-size distribution in 
these shells. They found from their sample of RSGs that 
R'-^ = 4.1 + 0.1E{B -V)~ 0.21og5, where log^ is the 
surface gravity of the RSG and has values of between 
—0.5 and 0.5 (when g is expressed in cgs units). 

We compared the observed SED for the progenitor 
to a model SED synthesized using STSDAS/SYNPHOT 
within IRAF from the MARCS mo del stellar atmo- 
spher es for RSGs at solar metallicity (jGustafsson et al.l 
I2008D with Teff in the range 3400-3800 K (essentially, 
spect ral types of M5 through late K; iLevesaue et al.l 
120051 ). in steps of 100 K, and surface gravities logg = 
—0.5, 0.0, and -1-0.5. The models are for 15 M0 stars, 
assuming spherical geometry and a microturbulence ve- 
locity of 5 km s~^. The filter response func t ion at J in 
the LCO system is taken from lPersson et all ()1998l ). and 
the function for Ks in the MKO system is from the on- 
line INGRID instrument pag£0. We allowed R'y to vary 
from the typical value of 3.1 to larger values and com- 
puted a r ange of Ay from 1.4 to 4.2 r nag fo r each R'y, 
following ICardelh. Clavton. fc Mathlsl (|1989D . We were 
able to eliminate the logg = -1-0.5 models outright, since 
none of these were able to reproduce the observed SED 
at any value of R'y, Ay, and Tcff- For all of the models 
at the other two surface gravities, we found that only the 
Teff = 3600 K models were allowed within the uncertain- 
ties in the photometry. Furthermore, we found that for 
the log g — —0.5 models, a range in R'y from 4.0 to 4.7 
was allowed, although the range in Ay was constrained 
to 3.05-3.10 mag. For the logg = 0.0 models, R'y could 
range from 4.1 to 4.6, and again, the extinction was con- 
strained to Ay = 3.10-3.15 mag. 

The total Ay in this case can be considered an excess 
extinction, due to circumstellar dust, together with the 
line-of-sight extinction to the star, which we assume to 
be the total extinction to the SN (see § 13. 5p . Since we 
have shown that the star was variable, possibly due to 
variations in the amount of circumstellar dust, this total 
Ay is essentially a time-averaged value over the duration 
of both the HST and ground-based observations. 

The average values of R'y and Ay (in mag) for the 
models at both surface gravities are 4.35 and 3.10, re- 
spectively. We show in Figure [3] a representative model 
at Toff = 3600 K with these values of R'y and Ay. The 
tightest constraint on the model SEDs comes from the 
/-band measurement; the large uncertain ty at , on the 
other hand, is not as constraining. From ILevesaue et al.l 

® http: / /www. ing.iac.es/Astronomy/instruments/ingrid/ingrid_filtcrs 



(|2005f) . this effective temperature would correspond to 

spectral type M3.^ 

Foll owing ILevesaue et al.l (j2005D and 

Bcssc U. Castelh. fc Ple^ (jl998f ). we computed the 
bolometric corrections at V, I, and from the 
MARCS RSG stellar atmosphere models at T^s = 3600 
K and at logg = —0.5 and logg = 0.0. We found 
BCy = -1.79, BCi = 0.33, and BC'k = 2.79 mag 
for logg = -0.5, and BCy = -1.78, BCi = 0.32, 
and BCks = 2.80 mag for logg — 0.0. Assuming 
that R'y = 4.35 and that Ay — 3.10 mag, we find that 
Mi,oi ^ -8.29±0.12, -8.29±0.11, and -8.13±0.22 mag 
from V, I, and Ks, respectively, adding in quadrature 
the uncertainties in the photometric measurements, in 
the inferred extinction (0.05 mag), in the bolometric 
corrections (0.01 mag), and in the host-galaxy distance 
modulus. It is very satisfying that our estimates of 
Mboi are exactly the same at V and /, which gives us 
confidence that the values for R'y, Ay, T^b, and logg, 
taken together, are all consistent. It also allows us 
to neglect Mboi obtained from the far less certain Kg 
measurement, although the value of Mboi from this band 
certainly agrees with those from the other two bands, 
to within the uncertainties. We therefore adopt the 
uncertainty-weighted mean, A/boi = —8.29 ± 0.08 mag, 
from V and /. We conservatively adopt an uncertainty 
of ±100 K in Toff, although the star's observed SED 
appears to tightly constrain the effective temperature of 
the star to less than this uncertainty. Assuming Mboi(0) 
= 4.74 mag, this corresponds to a bolometric luminosity 
relative to the Sun of log(Lboi/T0) = 5.21 ± 0.03. 

In Figure U] we show a Hertzsprung-Russell diagram 
with the locus of the SN 2012aw progenitor. For com- 
parison we also illustrate the massive-star evolutionary 
tracks at solar metallicity from Ekstrom et al. (2012|) for 
stars with initial rotation which is 40% of the critical ro- 
tation, at initial masses Mini = 15 and 20 M©. The star 
is clearly more luminous than the RSG terminus of the 
Mini = 15 M0 model, although the star's effective tem- 
perature is consistent with the value for the terminus of 
this model. The 20 M© model, however, terminates at 
a higher luminosity and far hotter Toff. The "red loop" 
of that track does approach the star's locus; however, 
based on the behavior of that track, we do not expect a 
star with this initial mass prematurely to reach its end- 
point at this luminosity and cooler effective temperature 
along the loop. We can infer, therefore, that the progen- 
itor's initial mass was in the range 15 < Mini(MQ) < 20. 
lEkstrom et al.l (j2012f ) do not provide model tracks be- 
tween Mini = 15 and 20 M©, so it is unclear at what 
luminosity and effective temperature stars with initial 
masses within this range would reach their termini. How- 
ever, by eye from the figure, we can interpolate that the 
star's locus could well be consistent with the endpoints 
of a putative 17 or 18 Mq track. Clearly, we require the 
actual model track to be more precise about this initial 
mass assignment. 

At the adopted effective temperature and luminosity, 
the star had an effective radiu s R = 1040 ± 100 Rq. 
From the lEkstrom et"al1 (|2012f ) Mi„i = 15 M© evolu- 
tionary track, the final mass is Mg„ = 11.1 M©, so the 
surface gravity would then be logg « —0.5. From the 
adopted Ay and R'y, the reddening E{B - V) = 0.71 
We then would expect, from the relation from 
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iMassev et all (|2005D between R'y, E{B - V), and \ogg, 
that R'y « 4.27. This indicates that there is consistency 
among these three parameters taken together, and this 
cross-check further provides us with confidence in our 
estimates of the star's properties. 

We note that the main difference between the anal- 
ysis we haye prese nted here and that presented by 
iFraser et al.l ()2012b[ ) — namely, that we are able to bet- 
ter constrain the effective temperature, luminosity, and, 
therefore, the initial mass estimate for the progenitor — 
stems not only from the differences in the photometry be- 
tween the two studies (we found that the star is brighter 
in V, and we were also able to reduce the uncertain- 
ties at J and Ks), but also from our assumption that 
Rv for the RSG progenitor was different from the typ- 
ical value of 3.1 for the diffuse ISM. One item to also 
note is that the evolutionary tracks (with out rotation; 
lEldridge. Izzard. fc Toul2008[) employed bv lFraser etlll 
( 2012bD tend to termin ate at significantly co oler effective 
temperatures than the lEkstrom et al.l ()2012|) tracks that 
we have used. The key, ultimately, was connecting the 
inferred properties of the SN 2012aw progenitor, based 
on our measurements, to tho se of the Galactic R SGs of 
similar luminosity and mass (|Massev et al.ll2005D . 

3.4. The Dust Around the Progenitor 

As'Mas sev et al.l (|2005( ) discussed, we can estimate the 
mass of the dust and duration of dust production that 
is responsible for the exc ess extinction which we infer 
for the RSG progenitor. IMassev et al.l assume a thin- 
shell approximation for the excess extinction, AAy = 
AR{3.2 X 103)Md/(47ri?2Ai?), where R is the stellar ra- 
dius (in m), AR is the extent of the thin dust layer 
or shell (or the path length through the dust) above 
the stellar surface, and is the dust mass (in kg). 
The circumstellar matter presumably extended several 
stellar radii above the star's surface, as in the case of 
a Ori, and the overall mass loss could well have been 
driven by convection in th e envelope (jJosselin fc Pled 
I2007t iChiavassa eraIll2010D . It is in the last ~ 1400 yr 
of the RSG phase, as seen in the 15 M© model from 
lEkstrom et al.l (|2012l ). when the star's luminosity, radius, 
a nd total mass loss pa rticularly increase. From Figure 4 
of IMassev etHI (1200511 . we see that the inferred AAy ~ 3 
mag (we will show in § 13.51 that the interstellar extinc- 
tion, for Rv = 3.1, is likely only Ay « 0.2 mag) is en- 
tirely consistent with the Mboi of the star. The dust 
production rate, M^, corresponding to this luminosity is 
~ 10~®'^^ Mq yr~^. From the relation above, the value 
of Md is ~ 6.1 X 10^^ kg. Comparing this to M^At, we 
find that a dust-producing episode of At « 1 yr, at some 
point prior to 1994 (since this dust already existed by the 
time of the first JfST/WFPC2 images), could account for 
AAy. Although quite a short interval of time, it is con- 
sistent with tha t inferred for episodic dust production in 
Galactic RSGs (jPanchi et al.lll99l iBester et al.l[T996l ). 

3.5. The Dust Around the Supernova 

However brief was the dust production, as lFraser et al.l 
(|2012b ) point out, the dust was far more quickly de- 
stroyed, likely by the X-ray / UV flash within hours of core 
collapse. As also noted by iFraser et al.l ()2012b[ ). such 
circumstellar dust destruction is not unprecedented for 



o ther SNe (e.g..lDwek et all [20081: [ Wesson et al.| [20Tol). 

iPoznanski et al.[ "( 20111 ) have found that the Na i D 
feature strength in low-resolution optical SN spectra is 
a poor in dicator of the amount of extincti o n to th e SN. 
However, IPoznanski . Proch aska. fc BloomI (|201 2b^ have 
established a well-calibrated relation between reddening 
and the equivalent width (EW) of the Na i D doublet, Dl 
and D2, based on more than one hundred high- resolution 
spectra of objects through a number of interstellar lines- 
of-sight. We have therefore measured the EW of the dou- 
blet, from both the Milky Way and host- galaxy compo- 
nents, as clearly detected in a high-resolution spectrum of 
SN 2012aw we obtained with th e High Resolution Echelle 
Spectrometer (jVogt et al.lll994D on the 10 m Keck I tele- 
scope in the red optics configuration ( "HIRESr" ) on 2012 
April 10.29. We used the C2 decker (i.e., the 0"86 slit), 
providing coverage of 3800-7300 A with a resolution 
of 50,000. The portion of this spectrum, centered around 
the Na i D feature, is shown in Figure [S] We find that for 
the Milky Way component, EW(D2 A5891.41) = 94 ± 8 
mA and EW(D1 A5897.39) = 56 ± 9 mA. For the com- 
ponent internal to M95, EW(D2 A5909.33) = 269 ± 14 
mA and EW(D1 A5915.32) = 231 ±11 mA. Uncertainties 
in the centroids of the absorption features are typically 
^ 0.05 A. (These are all vacuum wavelengths, corrected 
to a heliocentric frame of reference.) The features due 
to M95, particularly, the D2 component and, to a lesser 
extent, the Dl component, may be slightly saturate d. 

From the IPoznanski. Prochaska. fc BloomI (|2012bl ) re- 
lations, E{B - V) = 0.022 ± 0.013 mag from the 
Milky Way and E{B - T^) = 0.055 ± 0.014 mag for 
the host galaxy. (The uncertainties here are primar- 
ily from the systematic uncertainties in the relations; 
the measurement uncertainties are comparatively neg- 
ligible.) These relations take into account saturation in 
the features, which is a smaller effect compared to the 
systematic uncertainties in the relations. Our redden- 
ing estimate is consistent with the Galactic foreground 
redden ing estimate from iSchlegel. Finkbeiner. fc DavisI 
(fT998l ). E{B - V) = 0.028 mag. The total redden- 
ing indicated from the high-resolution SN spectrum, as- 
suming Rv = 3.1, is then E{B - V) = 0.077 mag, 
which is relatively low and comparable to the esti mated 
EjB - V) = 0.1 mag fouii d for SN 1999em (Baron^eial] 
[2OOOI; [Leonard et al.l [20021: lElmhamdi et al.. .2003). The 
visual ex tinction to the SN is then Av = 0.24 mag. 
Although IPoznanski et all (|2009l ) found that for a sam- 
ple of SNe II- P the best-fit i?,^ ~ 1-7, the relations 
from FPoznansk i. Prochaska. fc BloomI likely do not de- 
pend strongly on Rv, and therefore the reddening and 
extinction values we have estimated likely will not be 
significantly different if Rv 7^ 3.1, especially local to 
the SN. Nonetheless, it is evident that the SN explosion 
must have destroyed much or all of the circumstellar dust 
around the progenitor, leaving only what is most likely 
interstellar line-of-sight extinction, which we measured 
from the high-resolution SN spectrum. 

4. CONCLUSIONS 

The probable progenitor of the intermediate- 
luminosity SN II-P 2012aw/PTF12bvh has been 
identified in archival HST optical and ground-based 
near-infrared images. Using the photometry extracted 
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from those images we have constructed a SED for the 
star, and we analyze the SED to show that the star was 
a luminous (Mboi = —8.29 mag) RSG with spectral type 
M3 (Tcfi = 3600 K) and with substantial circumstellar 
dust up to 18 yr before explosion. This dust was 
clearly destroyed by the explosion, since the current 
extinction to the SN is relatively low. Although the 
existing, state-of-the-art, theoretical stellar evolutionary 
tracks do not terminate at the locus of the star in the 
Hertzsprung-Russell diagram, we surmise that the star 
had an initial mass ~ 17-18 M0. 

This mass is essentially the same as th e upper limit 
to the initial mass, 16.5 ± 1.5 Mq, that iSmartt et all 
(|2009D derived for SN II-P progenitors. It is not evident 
whether this result for SN 2012aw can be generalized 
for other intermediate-luminosity SNe II-P. The initial 
mass for the p rogenitor of SN 2004et is still a su bject 
of debate (e.g.. lLi et aT1l2005t [Crockett et al.ll2011( ). Ad- 
ditionally, the (7cr) upper lim it to the detection of SN 
1999em at / = 22.0 mag from ISmartt et all H2002f ). at a 
distance of 11.7 Mpc (jLeonard et al.ll2003l which implies 
a luminosity limit only a factor ~ 1.4 fainter than if SN 
1999em were at 10.0 Mpc, as is SN 2012aw), precludes 
detection of an analog to the SN 2012aw progenitor. It 
therefore remains uncertain what is the upper limit on 
the initial mass of the RSGs that give rise to "normal," 
intermediate-luminosity SNe II-P. The detected, unusual 
progenitors of high-lum inositv SNe II-P (e.g., SN 2008cn; 
lElias-Rosa et al.ll2009( ) may provide some indication of 



this limit. The rece nt stellar evolutionary tracks from 
lEkstrom et al.l ([20121 ), as well as those wi t h pul sationally 
driven superwinds bv lYoon fc CantieIlo[ ([2O1O0 . demon- 
strate the role of rotation and mass loss in the late-stage 
evolution of stars with Mini ^ 20 M0 . Still to be investi- 
gated more fully are the influences of factors such as the 



metallicity and binarity (e.g.. lSmith et al]|2011f) . To ver- 
ify the candidate SN 2012aw progenitor, we will need to 
return at very late times, when the SN has substantially 
faded, to see whether the dusty RSG has vanished. 
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TABLE 1 
Photometry of the SN 
2012AW Progenitor^ 



Band 


Magnitude 


F336W 
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F439W 


> 25.9'' 


F555W 


26.49(07) 


V 
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F814W 


23.39(02) 


/ 


23.44(02) 


J 


21.02(03) 


Ks 


19.47(19) 



Uncertainties (Ifr) are given 
in parentheses, in units of 
0.01 mag. 

'° These are 3 cr upper limits. 
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Fig. 1.— (a) A portion of the archival HST F814W image of M95 from 1994/1995; the star detected at the precise location of SN 2012aw 
is indicated by tick marks. The SN site is very near to the edge of the image mosaic, (b) A portion of the Kp-band AO image obtained 
using NIRC2 on the Keck-II telescope on 2012 March 27. 
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Fig. 2. — Photometry (mag) of the SN 2012aw progenitor measured with Dolphot from the individual HST/WFPC2 observations from 
program GO-5397 in the (a) F555W and (b) P^814W bands. Observations obtained on the same date in a given band have been averaged, 
weighted by the uncertainties in the individual measurements. Panel (c) shows the F555W— F814W color (mag). The dashed lines in all 
three panels arise from the uncertainty-weighted means of all measurements in each band, as returned by Dolphot. 
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Fig. 3. — The observed SED for the probable SN 2012aw progenitor. See Table [T] Also shown fo r comparison is an example SED 
computed via synthetic photometry of a MARCS model RSG stellar atmosphere (IGustafsson et ani2008l) with surface gravity logg = —0.5, 
effective temperature T^g = 3600 K, and extinction Ay = 3.10 mag (solid line), assuming an effective total-to-selective ratio of absorption, 
R'y = 4.35. Although a surface gravity logg = 0.0 and a limited range in R'^r are also allowed, both Ay and Toff are constrained by the 
observations; see the text. The model SED has been normalized at /, the band which provides the tightest photometric constraint. 
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Fig. 4. — Hertzsprung-Russell diagram showing the locus of the SN 2012aw progenitor. For comparison we also indicate massive-star 
evol utionary tracks at sola r metallicity, which include initial rotation on the main sequence that is 0.4 times the critical rotation velocity, 
from lEkstriim et al.l ||20T^ ') for Mini = 15 and 20 Mq. 



14 Van Dyk et al. 





n 


■"iiiivy Way 


M95 



5890 5895 5900 5905 5910 5915 



Heliocentric Wavelength (A) 



Fig. 5. — Spectrum of SN 2012aw obtained on 2012 April 10.29 with HIRESr on the 10 in Keck I telescope, centered on the Na I D 
absorption feature. Both feature components Dl A5896 A and D2 A5890 A are clearly detected from both the Galactic foreground (labeled 
as "Milky Way") and the SN host galaxy (labeled as "M95"). The features due to the host may be partially saturated. 



